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We present an optimal strategy having finite outcomes for estimating a single parameter of the 
displacement operator on an arbitrary finite dimensional system using a finite number of identical 
samples. Assuming the uniform a priori distribution for the displacement parameter, an optimal 
strategy can be constructed by making the square root measurement based on uniformly distributed 
sample points. This type of measurement automatically ensures the global maximality of the figure 
of merit, that is, the so called average score or fidelity. Quantum circuit implementations for the 
optimal strategies are provided in the case of a two dimensional system. 
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I. INTRODUCTION 

Quantum state estimation deals with how to estimate unknown parameters of a quantum state as precisely as 
possible. This problem was studied extensively in the 1970's in the context of the formalism of the probability 
: operator measure (POM) and of quantum Bayesian inference. Basic formulations were already established and 
' optimal strategies were found in various cases. The relevant works are reviewed in Refs. |Q,0|. 

In the 1990's, the same problem was revisited in a new context where one is allowed to use finite samples of a quantum 
system to be estimated, while most of earlier works were concerned with the estimation using a single sample. Massar 
and Popescu || obtained the optimal strategy for estimating a quantum pure state of a spin-1/2 system, say, \p), from 
N identically prepared samples \p)® N . Their strategy is based on the use of an infinite continuous set of projectors in 
the Hilbert space of \pf N . For a given unknown quantum system, the optimal estimation strategy is not unique. As 
O |' shown in Ref. one can always find an optimal strategy consisting of an infinite continuous set of POMs. On the 
other hand, for a finite dimensional system there must exist a discrete and finite POM achieving the same optimal 
bound as shown by Derka et al. 0]. From the point of view of the physical realization the latter is preferable, while 
the former might be a more useful mathematical tool to derive the maximum attainable average fidelity (which is 
^ a commonly used figure of merit for state estimation). In Ref. Q an algorithm is described for constructing such 
optimal and finite POMs for an arbitrary finite dimensional system in a pure state. Latorre et al. || then studied the 
optimal strategy with a minimum number of outcomes for a spin-1/2 system, and showed explicit forms of optimal 
minimal measurements for N = 1 ~ 5. Their analysis was extended to the cases of a mixed state of a spin-1/2 system 
||, and of an arbitrary spin system in a pure state JtJ. In Ref. ||, the closed form expressions for the maximum 
average score, the optimal minimal POM and its number of outcomes were derived by using the symmetric properties 
of the totally symmetric subspace supported by a tensor product of N identical samples. To construct the optimal 
minimal strategy explicitly, however, some parameters are to be determined and remain unsolved for larger N(> 6). 
For higher dimensional systems it becomes more difficult to find concrete forms for the optimal minimal strategy [[7|. 
As for the maximum average fidelity the explicit expression for an arbitrary N was obtained in Ref. Q . This bound 
was derived by using the fact that quantum optimal state estimation using N samples can be viewed as the limiting 
case M — > oo of universal optimal cloning generating M copies from N inputs for which the maximum average fidelity 
was given by Werner j|. [] 

Although a substantial progress has been obtained in the quantum state estimation, it is still a difficult and open 
problem how to find explicit and physically realizable solutions for optimal strategies in an algorithmic way, especially 
in the case of higher dimensional systems and larger numbers of samples. Moreover, discussions given so far in the 
literatures for ensuring the optimality of discrete and finite POMs were focused only on the condition for extremality 



1 The optimal cloning map of Ref. || may have a connection with the infinite continuous version of the optimal state estimation 
strategy. 
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and not on the full conditions for the existence of a global maximum, which are reviewed in Refs. |l],[2|. In general, 
seeking all extrema and picking up the point corresponding to the global maximum is not necessarily a trivial task 
for complex systems. 

In this paper, we focus on a single parameter estimation of an arbitrary finite dimensional system in a pure state 
and give finite element optimal strategies that can be constructed in a straightforward way and ensure the global 
maximality conditions for the POM. 



II. OPTIMAL PHASE ESTIMATION 



Consider a finite dimensional system described in a Hilbert space H and let {|0) , |1) , \ K)} be its basis built from 
the eigenstates of the observable O on 7i, O \k) = k \k). Such a system may be, e.g., an optical field produced by 
quantum scissors |Tof (with O the photon number operator) or a spin j system described by a superposition of the 
eigenstates {\j, m) , m = —j, j} (with O the spin operator). The problem we consider is the estimation of a unitary 

evolution u specified by a displacement parameter 6, that is, u(9) = e~ t0 ° . We suppose that the initial state of the 
system is a priori known and reads 

W0)>=5>l*> (2.1) 

fc=0 

(with cj; non zero arbitrary complex coefficients), but we do not have any a priori knowledge about 0. After the 
evolution, the system will be in a state 

K 

|^))=6(e)^(0))=^c fc e- iefc |fc). (2.2) 

fe=0 

It is assumed that N identical samples of the system are available. The combined system is then described on the 
totally symmetric bosonic subspace of 'H® N |^,^[ as 

= \mf N = Y^C{n)e- w T, k *»* \H) , (2.3) 

n 

where 

K nk 

C{n)^^N~\\{^= (2-4) 



fc=0 v 11 

means the summation over (K + l)-tuples ft = (no, ...,tik) with X^fcLo n k — N, and \n) (= \no, ...,iik)) is the 

( N + K\ 

occupation number basis. The dimensionality Db of this space is Db = ( ^ ]■ For our present purpose, however, 

it is enough to consider the smaller subspace spanned by the eigenstates of the compound operator Ot = 2Zi=i 0(i), 
where 0(i) is the observable for the i-th sample. Let {nf^} be the set of (K + l)-tuples that satisfy X)fc=o ^ Uk = J 



and define Aj = yY^i )l 2 - Then the state in Eq. ( |2.3| ) can be rewritten as 



KN 



\y(8))=J2Are- l9J \J), (2.5) 



./=n 

where 

|J>=A7 1 ^C(n,( J) )|n l (J) ), 6 T \J) = J\J). (2.6) 

i 

The POM describing the optimal estimation strategy is constructed in the Z?t = KN + 1 dimensional subspace Ht 
spanned by the set {| J)}. 

Such a POM {p, m } should maximize the following score 
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S(N) = J2^T I d6Tr((x m y(6))\(i> m \ip(e))\ 2 , (2-7) 

where $(6) = \ 1 i'{9)) (^(8) | and \ip m ) is a guessed state according to the m-th outcome of the measurement. Optimality 
can be discussed along with the conditions for quantum Bayesian optimization p],0]. Let us introduce the score 
operators 



W m = — d6^{6)\{^ m \^{6))\\ (2.8) 
27T Jo 

They include all the a priori information. Then the necessary and sufficient conditions such that a certain set {/2 m } 
globally maximizes the average score for a fixed set of {|VVn)} are expressed as 

(i) r = J2 m Wmfam is hermitian and (f - W m )jl m = 0, Vm, 

(2.9) 

(ii) t-W m >0, Vm, 

where T is called the Lagrange operator, and the average score (2.7) can then be rewritten as 

S(N) = Trf. (2.10) 

The optimal estimation strategy can be constructed in the following way. First take M states corresponding to 
uniformly distributed sample points in 9 € [0, 2w), that is, 

K 



\*l> m ) = ^c k <r i *& k \k) i (m = 0,.,M-l), (2.11) 



k=Q 



and let us denote its N tensor product states as |^ m ) = \ip m )® N ■ Then consider the state vector 



*-*|* m ), tf = £|tf m )(* m |. (2.12) 



The set {/t m = \fJ, m ) (£t m |} is easily seen to be a set of non-negative Hermitian operators satisf ying the resolution 
of the identity on TCt, and thus a POM. This POM is often called the square root measurement |pl]-[l3|. If we take 
M > KN + 1 sample points, {Ami also works as the optimal estimation strategy. Under the condition M > KN + 1 
we have, in fact, that (cf. Eq. (|2.5|) with 9 = 2-Km/M) 



KN 



# = Af^A 2 7 |J)(J|, (2.13) 

because 

M-l 

e i22 M Ln = MS nt o, for - KN <n< KN. (2.14) 



m=0 



Therefore, from Eq. ( |2.12| ) we get 



_L Ve-^ J |J). (2.15) 
To prove optimality, let us first rewrite the score function as (see Eqs. (^^) and (|2.1l| )) 

\(^ m \m)\ 2 =d + J2d L (e^-^ + e- 1 ^-^) , (2.16) 

L=l 

where 
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K-L 



d L= Y \ck+LCk\ 



k=0 



By substituting Eq. (2.16) into Eq. (|2.8|), we obtain 

KN K KN-L 

d Y A j\J)(J\ + Y dL Y AjA J+L (e^ L \J)(J + L\+e-^ L \J + L)(J\ 



We then have 



,7=0 



L=l 



,7=0 



M-l 



KN 



K KN-L 



t=J2w m fi m = d Q Y / A 2 j\J)(J\ + J2 d L Y A.jA J+L (\J)(J\ + \J + L)(J + L 



,7=0 



L=l 



,7=0 



(2.17) 



(2.18) 



(2.19) 



where the orthogonality relation Eq. ( 2.14 ) was used under the condition M > KN + 1 . The operator V — W m is now 

K KN-L 



L=l 



.7=0 



£-W m = Y,d L Y A J A W \J){J\ + \J + L){J + L\-(c^ L \J){J + L\+c-^ L \J + L){J\) 



(2.20) 



Now look at each operator enclosed by [...] in Eq. ( |2.20 ). This 2x2 matrix has the eigenvalues and 2 and is 
non- negative. So is the sum of them. Since the coefficients and Aj in Eq. ( 2.2Cj ) are also positive, the second 
condition (ii) of Eq. (2.9) is proved. The first condition (i) can be easily checked by direct calculation. The maximum 
average score is then given by 



A' A' 



K KN-L 



(2.21) 



7=0 



L = l 



7=0 



This maximum is independent of the number of sample points M, that is, for M > KN + 1 the attainable average 
score is saturated. In fact, in the limit of M — > oo we can construct the infinite continuous POM 



dH(ip) = \ii{tp)){n(<p)\dip/2ir, 



dn{ip) 



(2.22) 



where 



KN 

£< 

J=Q 



-iipj 



\J) 



(2.23) 



and t his at tains the same maximum as Eq. (2.21). Thus it is proved that the measurement state vector Eq. (2.12) 
(Eq. (2.15)) provides the optimal estimation strategy. For the minimum number of sample points M = KN + 1, 
{I Mm)} is an orthonormal set, that is, a von Neumann measurement. 

H ere we mention other strategies {A m } which ex trem ize the average score Eq. (2.7) (for the same W m as in Eq. 
( |2.8| )), that is, satisfy the first condition (i) of Eq. (2.9), but not necessarily the second condition (ii). Consider, for 
example, the states orthogonal to the vV-tensor-product sample states | \& m ) , that is, the states 



KN 



\<) = Y 



7=0 



KN 
J 



V(-l) J e 



| J). 



(2.24) 



The square root measurement for discriminating {^m)} i s given as 



M-l 



>=(El**X*m| 



m— 



| m/ 



KN 

= -^Y(-l) J e- i2 ^ J \J) 



(2.25) 



7=0 
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The Lagrange operator for this measurement is 

M-l KN K KN-L 

f^=Y / W m fii = doY / A 2 j\J)(J\ + Y, d t E AjAj +l (-1) l (\J)(J\ + \J + L)(J + L\), 

m=0 J=0 L=l J=0 

and we have that 



K KN-L 



(-l) L (\J)(J\ + \J + L)(J + L\) 



which is easily seen to satisly 



L = l j=o 
(e l2 ^ L \J){J + L\ + e~ l ^ L \J + L)(J\j 



(t 1 - - W m ) = 0. 



(2.26) 



(2.27) 



(2.28) 



Thus the POM {/i^ = (/x^J} is an extremal solution for the average score Eq. (2.7). However, e.g., in the case 
of a two dimensional system (K = 1 in Eq. (2.1)), Eq. (2.27) reduces to 



N-l 



f x - W m = -d! AjA J+1 



J=0 



\J)(J\ + I J + 1)(J + 1| + e l ^r\J){J + 1| + c- l ^\J+ 1) (J| 



(2.29) 



which is clearly a non-positive definite operator, that is, f x — W m < 0. In contrast to Eqs. ( |2.9| ), this means that 
{/t„} attains a global minimum of the average score. In the general case of K > 1, it is not necessarily the case that 
T 1 - — W m is positive or negative definite. Thus, in general, the set {ft^} ma y represent strategies that attain either 
local maxima or minima. 



III. QUANTUM CIRCUIT FOR THE OPTIMAL STRATEGY 



Let us consider physical implementations of the optimal estimation strategy represented by Eq. ( p. 15 ). This is 
a collective measurement on a finite sample system. This sort of measurement can, in principle, be realized by a 
quantum circuit acting on the combined system and a separable measurement on each sample system. From this 
point of view, the main problem is a synthesis of an appropriate quantum circuit. When we take the minimum 



number of outputs M = KN + 1, the measurement basis {|^t m )} is orthonormal, and Eq. (2.15) is the discrete Fourier 



transform in the subspace 7Yt- The discrete Fourier transform is a fundamental tool in quantum computation. The 
corresponding circuit working on qubit systems is already known (see, e.g., Ref. fl4|| ). Therefore in the case of a two 
dimensional system (that is, K = 1 in Eq. (2.2)), we may apply this result to synthesizing the optimal estimation 
strategy. 

Let us start with the simplest two dimensional system with N = 2. The state to be measured is 



\mf 2 = A |0) T + A lG ~ w |1) T + A 2 e~™ |2) 



T ' 



where Aq = Cg, A\ — v2cqCi, Ai = c\ , and 



|o) T =|oo), 

|1) T =-L(|01> + |10» : 
|2) T = 111) , 



(3.1) 

(3.2a) 
(3.2b) 
(3.2c) 



assuming that the coefficients c, are real for eliminating unessential phase factors. Here we have used the subscript T 
for denoting the basis | J) of the KN + 1 dimensional subspace Ht, and the basis states 1 00) , |01), 1 10) and |11) span 
the space TL® 2 . The optimal strategy with the minimal outputs is represented as 



|/x m ) = -±= (|0) T + |1) T + e'^ |2) T ) , (m = 0, 1, 2). 



(3.3) 
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This expression is, however, rather inconvenient for directly applying the discrete Fourier transform quantum network, 
which is usually defined in a 2 n dimensional space corresponding to n qubit systems. Therefore, it seems better to 
consider the measurement circuit in a 4 dimensional space H® 2 spanned by the basis { 1 00) , 1 01} , 1 10) , |11)}. In fact 
we also have the 4-output optimal, nonminimal (M = 4) strategy in 7ix as 

\fim) = -j= (|0) T +e- <2 ^ |l) T + e- 4 ^ |2) T ) , (m = 0,1, 2, 3), (3.4) 

and this can be extended to a von Neumann measurement in Ti® 2 . It is thus helpful to transform the state of the 
input samples by the unitary operator defined by the circuit shown in Fig. |l| such that 

t W \iP{6)) m = A |00) + A ie ~ ie |01) + A 2 e- ae |10) . (3.5) 

We may then apply the measurement 

| Mm ) = i(|00)+e- <2 ^|01)+e- i ^ i |10)), (m = 0,1,2,3), (3.6) 

or its Naimark extension in H® 2 

|n m ) = i=(|00)+e-^|01)+e-^|10) + e-^|ll)), (to = 0,1, 2, 3) (3.7) 

which, using the discrete Fourier transform J7dft shown by the circuit of Fig. ^, can be explicitly written as 

|n ) = UI FT |00) , (3.8a) 
iHr) = ^ FT |10), (3.8b) 
|H 2 ) = Ul FT |01) , (3.8c) 

|n 3 > =^ FT |ii). (3.8d) 



-(2) c 



H 



— © * 

FIG. 1. The circuit for the basis transformation from {|0) T , |1) T , |2) T } to { 00) , |01) , 1 10) } . H is the Hadamard transfor- 
mation. 



H 



DFT 



H 



FIG. 2. The circuit for the discrete Fourier transform on two qubit systems. The two bit gate R(<f>) performs the transfor- 
mation \x) \y) i-> e' 3 ^ \x) \y). 

Thus the optimal phase estimation can be realized by first performing the unitary transformation C/dftT^ 2 '' on the 
input state \i[>(9))® 2 , then measuring the transformed state in the basis {|00) , 1 10) , |01) , |11)} (which is a separable 
measurement), and finally deciding the phase as 6 = 0, f ,7r, or according to whether the outcome is |00), |10), 
1 01) , or 1 1 1 } , respectively. This is summarized in Fig. |^. 



•(2) 



DFT 



} 



00 
01 



1 1 





JL 
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10-*- n 

M 
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FIG. 3. The circuit structure for the optimal estimation strategy in the case of N = 2. 
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Similarly, in the case of TV = 3, the input state is represented as 

\^(9)f 3 = A |0) T + A 1 e~ ie |1) T + A 2 e~* 2e |2) T + A 3 e- m |3 

where 



T • 



|0) T - |000) , 

|1) T =i=(|001> + |010> + |100>), 

|2)x=^(|H0) 
|3) T = jlll) - 



1101) 



1011) 



(3.9) 

(3.10a) 
(3.10b) 

(3.10c) 
(3.10d) 



Let be the unitary operator which converts the basis states {|0) T , |1) T , |2) T , |3) T } into |Q)®{|00) , |01) , 1 10) , |11)}, 
respectively, that is 



f O) |^(#))® 3 = | ) ® (A |00) + A lC - i9 |01) + A 2 e~ l2e |10) + A 3 e- m |11)) . 



(3.11) 



The estimation strategy can then be co nstr ucted again in the 4 dimensional space 7i® 2 , where the minimal optimal 
measurement is actually given by Eq. (3.7) and can be realized just as in the previous case. The unitary operator 
T^ 3 ) can be effected by the circuit shown in Fig. |i] which consists of three main blocks. In the first block, the operator 
acts on the first two qubits of the Jir basis states, which gives 



|0) T 
|1>t 

|2> T 
|3) t 



|000) 

1 
i 

|101) 



1001) 
1 100) 



1010) 

ion) 



(3.12a) 
(3.12b) 

(3.12c) 
(3.12d) 



The second block transforms the state 1 100) into the state 1 110) and the last block S^ 3 \vx), which includes the 
conditional rotation on one qubit 



Vi 



l 



2 
3 

2 J_ 

3 V3 



(3.13) 



leads to the final basis states |0) <g> {|00) , 1 01) , 1 10) , 1 1 1) } as required. 

■ , - ; _' i 



H2) 



S (3) (v0 



S f3) (P,) 
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FIG. 4. The circuit which converts the basis states {|0) T , |1) T , |2) T , |3) T } into |0> ® {|00) , |01) , 10) , |11)}. 
Finally we mention the case of N = 4. Let us observe that the Ht basis can be written as 
0^ = |000)®|Q), 



|001) + |010) + |100) ) ® |0) + |000) ® |1) 



(3.14a) 
(3.14b) 
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5 > =4= 

1 



( |001) + |010> + |100> ) ® |i) + ( |no> + 1 101) + |011) ) ® |0) 



|110) + |101) + |011) ) <g> |1) + I 111) <g> |0) 



4j> = |111}<8>|1) 



(3.14c) 

(3.14d) 
(3.14e) 



The first 3 qubits can be transformed into |000), |001), |010), or |011) by applying f^K Thus the first qubit of all 
the basis states becomes |0), and can be factorized out. By further applying on the remaining 3 qubits a CC-NOT 
gate, the operator S^ 3 \ii2) with 



V2 



1 VI 

2 2 
\/3 1 

2 2 



(3.15) 



and the operator as shown in Fig. [|, the basis states { 

{|000) , |001) , |010) , |011) , |100)}, respectively. 



4) } arc finally transformed into 10) 
' T 



'(4) 



'(3) 



S (3 \v 2 ) 



W 



(3). 



w 



(3). 





T 

H 
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FIG. 5. The circuit which transforms the basis states {|o) , |4) } into |0)®{|000) , |001) , |010) , |011> , 1 100) }, respectively. 



Let T^ 4 ) be this basis transformation represented by the circuit of Fig. || After transforming the input state by 

it is then sufficient to perform the measurement on the last 3 qubits. Let {\L) 3 ; L = 0, 1, 7} be the 3 qubit basis 

{|000) , 1 001) , |111)}. Then the minimal optimal measurement is given as 



i 4 

_Y^ e -i^L\ L) ^ ( m = 0,1,2,3,4). 



(3.16) 



L=0 



However, for applying the discrete Fourier transform network, it is convenient to take the other optimal strategy 
consisting of the overcomplete states 



X 4 

J2e- i2 -^ L \L) 3 , (m = 0,l,...,7). 



L=0 



These can be orthogonalized in the 8 dimensional space 7i® 3 as 



i 7 

n ™) = 7f E e ~^ L l L ) 3 . ( m = o. !. 7 )- 



L=0 



DFT 



(3.17) 



(3.18) 



The circuit 



This is just the discrete Fourier transform in 7i® 3 , and can be written in the form |II„, 

corresponding to C/dft is found in Ref. |Q. Therefore the optimal estimation is realized by applying the unitary 
transform (/ eg) Udft) ■ on the input state |'0(0)) (84 and then by measuring the last 3 qubits of the transformed 
state in the basis {|L) 3 }. According to the outcome, we decide the phase to be exp(— i2irL/8), respectively. 

The case of larger N can be treated in a similar way by applying the circuits used in the case of lower N inductively. 
In the general case of higher dimensional systems, one should first develope basic tools for costructing quantum circuits 
(for some recent progress see, e.g., Ref. Q), and at present practical circuit synthesizations remain an open problem. 
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IV. CONCLUDING REMARKS 



We have shown how to construct the optimal strategies for estimating a displacement parameter of an arbitrary 
finite dimensional system in a pure state. These are based on the square root measurement for discriminating the 
states corresponding to the uniformly distributed sample points of the parameter. We have assumed that a priori 
probability distribution of the parameter is uniform. When the a priori distribution is not uniform, or the system to 
be estimated is in a mixed state, the strategy based on the square root measurement is not in general optimal. 

Within the assumption of a uniform a priori distribution and the purity of the system to be estimated, it is 
a remaining problem whether our method applies to the estimation of two or more parameters. The simplest case 
would be the estimation of a two state system using finite identical samples \ip(6, (f>))® N . As for the infinite continuous 
POM, we can show that the square root measurement 



d4>d9s'm9 
2^ ' 



o Jo 



da(e,4>) = i, 



(4.1) 



where 



KM)) 



2tt p-K 



d<t>' <ww (MM)) WM)lr 

4tt ./o Jo 



MM))®", 



(4.2) 



provides the optimal strategy. Whether its discrete and finite version should exist and be built from the uniformly 
distributed sample points of (6, </>) is still an open question at present. 
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